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Abstract 



Wc study the stochastic dynamics of growth and shrinkage of single 
actin filaments or microtubules taking into account insertion, removal, and 
ATP/GTP hydrolysis of subunits. The resulting phase diagram contains 
three different phases: a rapidly growing phase, an intermediate phase and 
a bound phase. We analyze all these phases, with an emphasis on the bound 
phase. We also discuss how hydrolysis affects force-velocity curves. The 
bound phase shows features of dynamic instability, which we characterize in 
terms of the time needed for the ATP/GTP cap to disappear as well as the 
time needed for the filament to reach a length of zero (i.e. to collapse) for 
the first time. We obtain exact expressions for all these quantities, which 
we test using Monte Carlo simulations. 

Key words: Actin; Microtubule; Dynamic instability; ATP hydrolysis; 
Force- velocity curves; Stochastic dynamics 
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Introduction 

A large number of structural elements of cells are made of fibers. Well 
studied examples of these fibers are microtubules and actin filaments. Mi- 
crotubules are able to undergo rapid dynamic transitions between growth 
(polymerization) and decay (depolymerization) in a process called dynamic 
instability yj). Actin filaments are able to undergo treadmilling-like mo- 
tion. These dynamic features of microtubules and actin filaments play an 
essential role in cellular biology (0). For instance, the treadmilling of actin 
filaments occurs in filopodia, lamellipodia, flagella and stereocilia (0, 0, S)- 
Actin growth dynamics is also important in acrosome reactions, where sperm 
fuses with egg (tS, iS, iS) . During cell division, the movements of chromosomes 
are coupled to the elongation and shortening of the microtubules to which 
they bind 0,0). 

Energy dissipation is critical for these dynamic non-equilibrium features 
of microtubules and actin. Energy is dissipated when ATP (respectively 
GTP) associated to actin monomers (respectively tubulin dimers) is irre- 
versibly hydrolyzed into ADP (respectively GDP). Since this hydrolysis 
process typically lags behind the assembly process, a c ap of consecutive 
ATP / GTP subunits can form at the end of the filament HO). 



Let us first consider studies of the dynamic instability of microtubules. 
Non-equilibrium properties of microtubules have often been described using 
a phenomenological two-state model developed by Dogterom and Leibler 
(jl2l ). In such a model, a microtubule exists either in a rescue phase (where 
a GTP cap exists at the end of the microtubule) or a catastrophe phase 
(with no GTP cap), with stochastic transitions between the two states. 
A limitation of such a model is that a switching frequency is built in the 
model rather than derived from a precise theoretical modeling of the GTP 
cap. This question was addressed later by Flyvbjerg et al. (jl3l . 114 ). where 
a theory for the dynamics of the GTP cap was included. At about the same 
time, a mathematical analysis of the Dogterom-Leibler model using Green 
functions formalism was carried out in ref. The study of Flyvbjerg et 



al. yj, 1 14 ). was generalized in (jld ). with the use of a variational method 



and numerical simulations. This kind of stochastic model for the dynamic 
instability of microtubules was further studied by Antal et al. (0,13). The 
Antal et al. model takes into account the addition and hydrolysis of GTP 
subunits, and the removal of GDP subunits. Exact calculations are carried 
out in some particular cases such as when the GDP detachment rate goes to 
zero or infinity, however no exact solution of the model is given for arbitrary 
attachment and detachment rates of both GTP and GDP subunits. 
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It was thought for a long time that only microtubules were able to un- 
dergo dynamic instability. Recent experiments on single actin filaments, 
however, have shown that an actin filament can also have dynamic-instability- 
like large length fluctuations ([l9j, 1201 ) . A similar behavior was observed in 
experiments where act in p olymerization was regulated by binding proteins 
such as ADF/Cofilin ([ill, l^i)- Vavylonis et al. ^ have studied theo- 
retically actin polymerization kinetics in the absence of binding proteins. 
Their model takes into account polymerization, depolymerization and ran- 
dom ATP hydrolysis. In their work, the ATP hydrolysis was separated into 
two steps: the formation of ADP-Pi- actin and the formation of ADP-actin 
by releasing the phosphate Pi. Vavylonis et al. have reported large fluctua- 
tions near the critical concentration, where the growth rate of the filament 
vanishes. More recently Stukalin et al. (0) have studied another model for 
actin polymerization, which takes into account ATP hydrolysis in a single 
step (neglecting the ADP-Pi-actin state) and occurring only at the inter- 
face between ATP-actin and ADP-actin (vectorial model) or at a random 
location (random model). This model too shows large fiuctuations near the 
critical concentration, despite the differences mentioned above. Note that 
both mechanisms (vectorial or random) are still considered since experiments 
are presently not able to resolve the cap structure of either microtubule or 
actin filaments. 

In this paper, we study the dynamics of a single filament, which can be 
either an actin filament or a microtubule, using simple rules for the chemical 
reactions occurring at each site of the filament. The advantage of such a 
simple coarse-grained non-equilibrium model is that it provides insights into 
the general phenomenon of self-assembly of linear fibers. Here, we follow 
the model for the growth of an actin filament developed in ref. (|24l). We de- 
scribe a new dynamical phase of this model, which we call the bound phase 
by analogy with two-state models of microtubules 

CI). 

The characteriza- 
tion of this bound phase is particularly important, because experimental 
observations of time-independent filament distribution of actin or micro- 
tubules (0) are likely to correspond to this phase. In addition, we analyze 
the dynamic instability with this model. We think that dynamic instability 
is not a specific feature of microtubules but could also be present in actin 
filaments. We argue that one reason why dynamic instability is less often ob- 
served with actin than with microtubules has to do with the physical values 
of some parameters which are less favorable for actin than for microtubules. 
This conclusion is also supported by the work of Hill in his theoretical study 
of actin polymerization (fiol . Ell ). In these references, a discrete site-based 
model for a single actin filament with the vectorial process of hydrolysis is 
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developed, which has many similarities with our model. 

In short, the model studied in this paper presents three dynamical phases, 
which are all non-equilibrium steady states: a bound phase (phase I), where 
the average cap length and the filament length remain constant with time, 
an intermediate phase (phase II) , where the average cap length remains con- 
stant and the filament grows linearly with time, and a rapidly growing phase 
(phase III) where the cap and the filament both grow linearly with time. 
The phases II and III were already present in the study of ref. (0), but 
phase I was not analyzed there. Thus the description of the main features 
of phase I (such as the average length, the distribution of lengths) is one of 
main results of this paper. 

In addition, we discuss how GTP/ATP hydrolysis affects force- velocity 
curves and we characterize the large fluctuations of the filament, by calcu- 
lating the time needed for the cap to disappear in phase I and II as well as 
the time needed for the complete filament to reach a length of zero {i.e. to 
collapse) for the first time in phase I. Due to the simplicity of the model, 
we are able to obtain exact expressions for all these quantities. We also test 
these results using Monte Carlo simulations. 



Model 

We study a model for the dynamics of of single actin or microtubule filaments 
taking into account ATP/GTP hydrolysis. Our model is very much in the 
same spirit as that of ref. (0) and has also several common features with the 
Hill et al. model 111). We assume that polymerization occurs, for actin. 



via the addition of single ATP subunit (GTP subunit for microtubule), at 
the barbed end (plus end for microtubule) (jil) of the filament. We assume 
that the other end is attached to a wall and no activity happens there. 
Let U and Wt be the rates of addition and removal of ATP / GTP subunits 
respectively, which can occur only at the filament end. The subunits on the 
filament can hydrolyze ATP/GTP and become ADP/GDP subunits with 
a rate R. We assume that this process can occur only at the interface of 
ATP-ADP or GTP-GDP subunits. This corresponds to the vectorial model 
of hydrolysis, which is used in the Hill et al. model (10, ll|). Once the whole 



filament is hydrolyzed, the ADP/GDP subunit at the end of the filament can 
disassociate with a rate Wd. The addition, removal and hydrolysis events 
are depicted in Fig. [TJ We denote by d the size of a subunit. 

This model provides a simple coarse-grained description of the non- 
equilibrium self-assembly of linear fibers. More sophisticated approaches are 
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possible, which could include in the case of actin, for instance, additional 
steps in the reaction such as the conversion of ATP into ADP-Pi-actin or 
the possibility of using more than one rate for the addition of ATP-subunits. 
It is also possible to extend our model to include growth from both ends of 
the filament rather than from a single end, as discussed in ref. (jj^ ). Another 
feature of actin or microtubule filaments, which we leave out in our model, is 
that these fibers are composed of several protofilaments (2 for actin and typ- 
ically 13 for microtubules). In the case of actin, it is reasonable to ignore the 
existence of the second protofilament due to strong inter-strand interactions 
between the two protofilaments (|24l ). In fact, we argue that the model with 
a single filament can be mapped to a related model with two protofilaments 
under conditions which are often met in practice. Indeed the mapping holds 
provided that the two protofilaments are strongly coupled, grow in paral- 
lel to each other and are initially displaced by half a monomer. The two 
models can then be mapped to each other provided that d is taken to be 
half the actin monomer size d = 5.4nm/2 = 2.7nm. This mapping suggests 
that many dynamical features of actin should already be present in a model 
which ignores the second protofilament. Similarly, microtubules may also 
be modeled using this simple one-filament model, in a coarse-grained way, 
provided d = 8nm/13 = 0.6nm is equal to the length of a tubulin monomer 
divided by 13, which is the average number of protofilaments in a micro- 
tubule (0). Keeping in mind the fact that the present model is applicable 
to both actin and microtubules, we use a terminology appropriate to actin 
to simplify the discussion in the rest of the paper. 

The actin filament dynamics is studied in terms of two variables n, the 
number of ADP subunits, and k, the number of ATP subunits, as shown 
in figure [TJ The dynamics of this system may be represented as a biased 
random walk in the upper quarter 2D plane (n, k). For instance, the addi- 
tion of one ATP subunit with rate U corresponds to a move in the upward 
direction. The removal of ATP subunits with Wt corresponds to a move 
in the downward direction. The hydrolysis of an ATP subunit results in an 
increase in n and decrease in k, both by one unit, which corresponds to a 
move in the diagonal direction as shown in the figure. The removal of ADP 
subunits can happen only when the cap is zero and therefore corresponds 
to a leftward move along the k = line. Let -P(n, k, t) be the probability 
of having n hydrolyzed ADP subunits and k unhydrolyzed ATP subunits at 
time t, such that I = (n + k)d is the total length of the filament. It obeys 
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Figure 1: (a) Schematic diagram showing addition with rate U, removal with 
rates Wt and Wd, and hydrolysis with rate R: (i) for the case where the 
cap length is non-zero and (ii) for the case where the cap length is zero. T 
stands for ATP (GTP) bound actin (microtubule) subunits while D stands 
for hydrolyzed ADP(GDP) subunits. At the T-D interface shown in (i), the 
hydrolysis occurs with a rate R. (b) Equivalent representation of the model 
as a biased random walk in the upper quarter 2D plane with rates U to go 
up (north), Wt to go down (south), R to go south-east and Wd to go west. 
The Wd move is only possible along the k = (southern boundary) line. 
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the following master equation: For A; > and n > we have 

^' = UP(n,k-l,t)+WTP(n,k+l,t)+RP(n-l,k+l,t)-(U+WT+R)P(n,k,t). 

at 

(1) 

When n = in Eq. [H P{—1, k + l,t) is set equal to zero. 
For k = and n > 1 we have, 

dP(n t) 

-^^^ = WoPin+l, 0, t)+WTP{n, 1, t)+RP{n-l, 1, t)-(U+WD)P{n, 0, t). 

(2) 

If /c = and n = 0, we have 

"^^^"'^"'^^ = WtP{0, l,t) + WdP{1, 0, t) - UP{0, 0, t). (3) 
The sum of the probabilities is normalized to 1 such that 

oo oo 

j;j;P(n,fe,t) = l. (4) 

n=Ofc=0 

We define the following generating functions 



Fk{x,t) = ^P{n,k,t)x'', (5) 

n>0 

Hn{y,t) = Y.P{n,k,t)y\ (6) 

fc>0 

G{x,y,t) = j;j;P(n,fe,t)xV (7) 

n>Ofc>0 

All quantities of interest can be computed from G{x,y,t). For instance, 
the average length of the filament is 

(0 = [(n> + Wl d = . {^^)^_^ + " {^^)^_^ • 

where (fc) is the average number of the ATP subunits and (n) is the average 
number of the ADP subunits. The average velocity of the filament is 

V = lim — (9) 

t— ♦oo at 

and the diffusion coefficient is 
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Similarly, the average velocity of the cap is 

d{k) 



J 



lim 



d, 



t—>-oo dt 

and the diffusion coefficient of the cap is 

D, = lim -4- {{k^) - {kf) 
t^oo 2 dt^^ ' ^ ' ' 



(11) 



(12) 



The solution of Eqs. [H [21 and [3] and the details of the calculations are given 
in the Appendix A. 



Phase diagram 

Our model for the dynamics of a single filament with ATP hydrolysis leads 
to the following steady-state phases: a bound phase (phase I), an interme- 
diate phase (phase II) and a rapidly growing phase (phase III). In phase I, 
the bound phase, the average velocity of the filament vi and the average 
velocity of the cap both vanish. Thus, the average filament and cap lengths 
remain constant in the long time limit. In phase /I, the filament is growing 
linearly in time, with a velocity vn-, but the average ATP cap length re- 
mains constant as a function of time. In phase //I, the filament as well the 
ATP cap are growing linearly in time with a filament velocity vm and cap 
velocity J. The boundary between phases I and II is the curve of equation 
vii = 0, and the boundary between phase II and III is the curve of equation 
J = 0. 

We have carried out simulations of the dynamics of the length of the 



filament, using the Gillespie algorithm (1251 ). According to this algorithm, 
the time to the next on-, off-, or hydrolysis-event is computed stochastically 
at each step of the simulation. We find that our simulation results agree 
with the exact calculations. 



The Bound phase (Phase I) 

In the representation of the model as a biased random walk shown in figure[Tl 
there is a regime of parameters for which the biased random walker converges 
towards the origin. After some transient time, the random walker enters 
what we call a bound steady state, where the motion of the walker is confined 
to a bounded region containing the origin. In the representation of the model 
as a filament, the filament length fluctuates as function of time around a 
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Figure 2: Filament length as a function of time for the three different phases 
of the model, as computed using Monte Carlo simulations. The value of 
the rates which were used are given in Table [TJ Figure (a) represents the 
bound phase (phase I), where a structure of avalanches in the evolution of 
the length can be seen. These avalanches correspond to series of sudden 
depolymerization events (collapse) followed by slow polymerization events 
(rescue). Figure (b) represents the intermediate phase (phase II) and the 
rapidly growing phase (phase III). In the intermediate phase (II), the dy- 
namics show large length fluctuations as compared to the rapidly growing 
phase (III), where the length fluctuations can hardly be resolved. 



time-independent average value (l) and at the same time, the cap length also 
fluctuates as function of time around a different time-independent average 
value {k)d. A typical evolution of the total length of the filament l{t), 
obtained from our Monte Carlo simulations, is shown in Fig. [2j 

We first discuss the properties of the cap before considering that of the 
total length. In the steady state {t oo), Fii.{x = 1) represents the dis- 
tribution of cap lengths, as defined in Eq. \5\ As shown in the Appendix 
A, 

Ffc(x = l) = (l-g)g^ (13) 



where 



U 



Since Fq = \ — q, -we see that q has the meaning of the probability of finding 
a non-zero cap in the steady state (|24| ). We consider for the moment only 
the case g < 1, which corresponds to phases I and II. From -Ffc(x = 1), we 
find that the average number of cap subunits is given by 



(15) 
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290 


4 
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R, ko and subunit length 



2J) and they characterize 



Table 1: Numerical estimates of the rates Wt, Wr 
d. The parameters for actin are taken from refs. (j^ 
the barbed end of an actin filament. The parameters for microtubule hold 
similarly for the plus end and at 35° C. The rate constant ko and rate Wo 
are taken from (ji), while the other two rates, Wt and R, are deduced from 
fitting the zero-force velocity data of ref. (fj^ l and the critical concentration 
measurements of ref. (j27l ). 



and0 



_Q_+ q 
(1 



^2- 



(16) 



As expected, these quantities diverge when approaching the transition to 
phase III when g — > 1. The standard deviation of the cap length is 



(7^ 



{k') - {kf 



(l-g)2- 

The relative fluctuations in the cap size are large since 



(17) 



{kf 



> 1. 



(18) 



We now investigate the overall length of the filament in the bound phase. 
This quantity together with the distribution of length in the bound phase 
can be obtained from the time-independent generating function G{x, y). The 
details of the calculation of this quantity are given in appendix A. We find 



G{x,y) 



(y - y+)[WT{y- - l)x + Rx{y^ - x) + WoV-il - x)] 
where y± are defined by 
1 



(19) 



y± = ^(u + WT + R± Vi.U + Wt + RY - AU{Wt + Rx)j . (20) 



^ Note that this expression of (k^) differs from that found in ref. (|24l). This discrepancy, 
we beheve, is probably due to a misprint in ref. |23 ). 
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From the derivatives of G{x,y), we obtain analytically the average length 
(/) using Eq. [8] as 



(0 



q {R^ + WdR^ + 2R^Wt + W^R + 2WdWtR + WdW^) 
Vii{q-1){WT + Ry 



vn{q - 1){Wt + R) 



(21) 



where vn is the shrinking velocity (since vn 
filament 

^\J{Wd^R) 



< in this regime) of the 



Vll 



Wd 



d. 



(22) 



Wt + R 

The length (/) diverges since vn when approaching the transition line 
between phases I and phase II. The length (/) as given by Eq.[2T]is plotted in 
Fig. [3] for the parameters of table I. We compare this exact expression with 
the result of our Monte Carlo simulations where the average is computed us- 
ing 1000 length values taken from different realizations. Excellent agreement 
is found with the analytical expression of Eq. [2TJ According to a simple di- 
mensional argument, the average length (/) should scale as —Djj/vfj, where 
Dij and vjj are the diffusion coefficient and velocity of phase II. We find 
that this scaling argument actually holds only close to the transition point 
between phase I and II. On the boundary line between phases I and II, the 
average filament velocity vanishes and hence the filament length is effectively 
undergoing an unbiased "random walk" . In such a case, we expect that on 
the boundary line {l'^) ~ t. We have also considered the fluctuations of l{t) 
using the standard deviation a defined as 



a 



{I') - {If 



(23) 



for which an explicit expression can be obtained from G(x, y). In Fig. [3l o" is 
shown as a function of U . Note that a is larger than (/), which corresponds 
to dynamic- instability- like large length fluctuations. 

In the limit ^ 0, ATP hydrolysis can be ignored in the assembly 
process. The model is then equivalent to a simple ID random walk with 
rates of growth U and decay Wt- In this case, phases II and III merge 
into a single growing phase. We find from Eq. [21] that {I) = Ud/iWr — U), 
which diverges as expected near the transition to the growing phase when 
U ~ Wt- According to the simple dimensional argument mentioned above, 
this length must scale as —Djv in terms of the diffusion coefficient and 
velocity of the growing phase (Il4l ). This is the case, since D = d^{U + Wt) /2 
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Figure 3: Average length (1) and its standard deviation a in the bound 
phase for (a) actin, and (b) microtubule as a function of the insertion rate, 
U (lower X-axis) , and of the concentration of the free ATP / GTP sub units C 
(upper X-axis). The solid lines represent analytical expressions of (/) (lower 
curves) and a (upper curves). Filled symbols represent values obtained from 
simulations for these quantities. All the curves and symbols are obtained 
using values of the rates given in Table. [TJ 



and V = (U — WT)d < and thus (/) = Ud/{WT — U) near the transition 
point. 

We have also computed the filament length distribution, P{1) in this 
phase using Monte Carlo simulations, as shown in Fig. HI In the inset, we 
compare the numerically obtained distribution with the following exponen- 
tial distribution 

P(/) = Poexp (-//(/)), (24) 

where (Z) is given by Eq. [2TJ In this figure, the distribution appears to be 
close to this exponential distribution. For any exponential distribution, the 
standard deviation, a should equal the mean {I). However as seen in figure 
[3l there is a difference between a and (/). Hence the distribution can not be 
a simple exponential, which could also have been guessed from the fact that 
the expression of G{x, y) is complicated. The exact analytical expression of 
the distribution could be calculated by performing an inverse Z-transform 
of the known G(x, y). 

In the bound phase, experiments with actin in ref. (figl ) report an average 
length in the 5-20 //m range at different monomer concentrations, and ex- 



periments with microtubules of ref. (|27l ) report a range 1-20 //m at different 



temperatures. Neither experiment corresponds precisely to the conditions 
for which the rates of table [1] are known. Thus a precise comparison is not 
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Figure 4: Filament length distribution in the bound phase for actin. (a) 
For U=1.5 and (b) for U=1.6 with all other parameters taken from 
Table [TJ The inset shows the same quantities {I vs P(l)) with P{1) in the 
log scale. The straight line in the inset is given by Eq. | 



possible at the moment, although we can certainly obtain with the present 
approach an average length in the range of microns using the rates of table 
[T] as shown in figure [3l 

Intermediate phase (Phase II) 

In the intermediate phase (phase II), the average ATP cap length remains 
constant as a function of time, while the filament grows linearly with time. 
The presence of this cap leads to interesting dynamics for the filament. A 
typical time evolution of the filament length is shown in figure [2j One can 
see the filament switching between growth (polymerization) and decay (de- 
polymerization) in a way which is completely analogous to what is observed 
in the microtubule dynamics (12). 

In this phase II, the average velocity of the filament is 

vn =[U -WTq-WD{l-q)]d. (25) 

This expression of vjj is the same as Eq. [22] except that now Eq. [25] corre- 
sponds to the regime where vjj > 0. The diffusion coefficient in this phase 
is 



Dii = — 
2 



U + WTq + WD{l-q) + 



2{Wd-Wt){U + WDq) 
Wt + R 



(26) 
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The expressions of vu and Djj are derived in the appendix A. When d equals 
to half the size of an actin subunit, we recover exactly the expressions of 
ref. (0). 

The transition between the bound phase (I) and the intermediate phase 
(II) is delimited by the vu = curve. When going from phase I to II, the 



average length in Eq. [21] varies as iU — Wd 



and the variance 



Wt+R ' ' 
-2 

. The transition from the 



of the length varies as {u - Wd 
intermediate phase II to the rapidly growing phase III is marked by a similar 
behavior. The cap length diverges as {U — Wt — R)~^ ■, and the variance of 
the fluctuations of the cap length diverges as {U — Wt — 



Rapidly Growing Phase (Phase III) 

In phase III, the length of the ATP cap and that of the filament are growing 
linearly with time. Thus the probability of finding a cap of zero length is 
zero in the limit t — > oo, that is Fq[x = 1) = 0. This also means that the 
probability of having a filament of zero length is also zero, i.e. P(0,0)=0. In 
this case, Eq. [49] reduces to 



dG{x,y,t) 
dt 



Using Eqs. 
following quantities 



C/(x-l)+I^rQ-l) G{x,y,t). (27) 

)5] and [SU] of Appendix A, one can easily obtain the 



VIII 

Dm 
J 



[U - Wt] d, 
d^ , 

y(^ + ^^T), 

[[/ - {Wt + R)\ d, 
d^ 

Y{U + Wt + R) 



(28) 

(29) 
(30) 

(31) 



Note that these quantities can be obtained from Eqs. [25][26l by taking 
the limit g — > 1, which marks the transition between phase III and phase II. 
In Fig. [2] the filament length is plotted as a function of time. Note that in 
this phase the velocity and the diffusion coefficient are the same as those of 
a filament with no ATP hydrolysis. The physical reason is that in phase III, 
the length of the non-hydrolyzed region (cap) is very large and the region 
with hydrolyzed subunits is never exposed. 
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Effect of force and actin concentration on active 
polymerization 

The driving force of self-assembly of the filament is the difference of chemical 
potential between bound and unbound ATP actin subunits. Since the chem- 
ical potential of unbound ATP actin subunits depend on the concentration 
C of free ATP actin subunits and on the external applied force /, the rates 
should depend also on these physical parameters. In the biological context, 
this external force corresponds to the common situation where a filament 
is pushing against a cell membrane. For the concentration dependance, we 
assume a simple first order kinetics for the binding of ATP actin monomers 
given that the solution is dilute in these monomers. This means that the 
rate U of binding of ATP actin is proportional to C while Wt, Wd, and R 
should be independent of C (jisl . [29I . Isol ). For the force dependance of the 
rates, general thermodynamical arguments only enforce a constraint on the 
ratio of the rates of binding to that of unbinding (0, Isil ) . A simple choice 
consistent with this and supported by microtubule experiments (3^) is to 
assume that only the binding rate i.e. U is force dependent. A more sophis- 
ticated modeling of the force dependance of the rates has been considered 
for instance for microtubules in (jsil ). All the constraints are then satisfied 
by assuming that U = koC exp {—fd/kBT), with ko, Wt, Wd, and R all 
independent of the force / and of the concentration C. We assume that 
/ > 0, so that the on-rate is reduced by the application of the force. The 
cap velocity in the rapidly growing phase(phase III), given by Eq. [30l can 
be written in terms of / and C as 

C) = \koCe^f'^/''B^ - {Wt + R)] d. (32) 

The phase boundary between phase II and phase III is defined by the curve 
J(C, f = fc) = 0. Equating the cap velocity to zero, we obtain the charac- 
teristic force, 



Wt + R 



ksT Co 

■^In-. (33) 



koC 

where the concentration Co is defined as 

Co = {Wt + R)/ko. (34) 



Below fc, the system is in phase III. This is also the point where q = 1. 
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The force-velocity relation in the intermediate phase is rewritten, using 
Eq. [251 as 



vnif,C) = koCe-f''/''^^ 



Wp + R 
Wt + R 



d - Wod. 



(35) 



The stall force fs is by definition the force at which vu{f = fs,C) = 0. 
From Eq. \35\ we obtain 



d 



In 



Wt + R\ Wd 



Wd + RJ koC 



(36) 



which can be written equivalently in terms of the critical concentration of 
the barbed end Ccrit as 



fs 



ksT ( Ccrit 



where 



crit 



Cn 



< Cn 



(37) 



(38) 



Wd + R, 

In the absence of hydrolysis, when i? ^ 0, we have Ccrit = Co and Eq. [371 
gwes the usual expression of the stall force given in the literature Q,, 23,, <2^, 

The velocity of the filament is shown in Fig. [5l with a zoomed in version 
marked as (b). This figure shows that for / < fc, the filament is in phase 
III, and that there the velocities of the filament with ATP hydrolysis or 
without are the same. At f = fc, the force- velocity curve changes its slope, 
as shown after the vertical line in figure [5l (see the zoomed- in figure [5)3) . 
When the concentration rather than the force is varied, a similar change of 
slope is observed at C = Co, which is accompanied by a discontinuity of the 
diffusion coefficient slightly above the critical concentration (j2l ). 

For fc < f < fs, the filament is in the intermediate phase, where the 
velocities in the presence and in the absence of ATP hydrolysis differ. The 
stall force with ATP hydrolysis is smaller than that in the absence of ATP 
hydrolysis. In view of this, a useful conclusion is that it is important to take 
into account the ATP hydrolysis for estimating the velocity of a filament 
when the force is close to the stall force. 

For f > fs, the velocity of the filament vanishes. It must be noted 
that, in this phase, the instantaneous velocity can be positive or negative, 
but the average velocity, in the long time limit, is zero. Another important 
point to note is that when the filament is stalled, ATP is still hydrolyzed. 
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. . Li I I 1 u u 1 1 

1 2 3 4 2.7 2.8 2.9 3 3.1 3.2 3.3 

Force (pN) Force (pN) 



Figure 5: (a) Steady-state force- velocity relation for a single actin filament 
shown for C = 1/iM (solid curve), (b) Zoom of the force- velocity relation 
near the stall force. The vertical lines represent fc and fs as shown in figure 
(b). For / < fc, the filament is in the rapidly growing phase (III) and the 
velocity is given by Eq. [28l In the bound phase (I) the velocity is zero. In 
the intermediate phase the velocity is given by Eq. [25l The dash-dotted line 
(figure (b), phase I and II) is given by Eq. [28l showing that the stall force 
is higher when ATP hydrolysis is neglected. 



This is analogous with models of molecular motors containing more than 



one cycle (|35|, l36|, 1371 ). Including the chemical cycle of ATP hydrolysis in 



addition to the mechanical cycle of addition/removal of subunits is for this 
reason important in the context of actin and microtubule models. One could 
imagine testing these predictions on the effect of ATP hydrolysis on force- 
velocity relations by carrying out force-velocity measurements near stalling 
conditions of abundant ATP or when ATP is sequestered by appropriate 
proteins (jsi ). 

All these observations can be summarized in a phase diagram in the 
coordinates / and C as shown in Fig. [6l As shown in Fig. [6l when C < Cq, 
the filament is either in the intermediate phase (II) or in the bound phase (I). 
In this region of the phase diagram, the fluctuations of the filament length 
are large as compared to the very small fluctuations observed in the rapidly 
growing phase (III). The large fluctuations observed in phase I correspond 
to the dynamic instability. 

In the case of microtubules, we find as shown in Tabled Cq = 8.75^M. 
This value is rather large when compared to typical experimental concentra- 
tions C ~ 1 — lO^M for microtubules. Thus microtubules are usually found 
in phases I and II where the length fluctuations are large and dynamic in- 
stability is commonly observed. 
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Figure 6: Phase diagram as a function of the normalized force fd/ksT and of 
the log of the ratio of the ATP subunit concentration C to the characteristic 
concentration Cq for actin. The phase diagram shows the bound phase 
(phase I), the intermediate phase (phase II), and the rapidly growing phase 
(phase III). The boundary hue between phase II and phase III is the curve 
J = and the boundary line between phase I and phase II is the curve 
VII = 0. 





Co ifxM) 


Ccrit (^M) 


fc (pN) 


fs (pN) 


Actin 


0.147 


0.141 


2.916 (at 1 ^M) 


2.978 (at 1 ^M) 


Microtubule 


8.75 


8.63 


5.65 (at 20 fiM) 


5.74 (at 20 ^M) 



Table 2: Estimates of the characteristic concentrations Cq and Ccrit at zero 
force, and of the characteristic forces fc and fs at a concentration of 1 /uM 
for actin and 20 /iM for microtubule using the rates of Table [H 



In the case of actin, we find that Cq = 0.147 /xM. Typical experimental 
actin concentrations are above this estimate, therefore, at zero force, actin 
filaments are usually seen in phase III. This may explain why the dynamic 
instability is rarely seen in actin experiments with pure actin. However, at 
small concentrations close to the critical concentration, large length fluctu- 
ations have been observed for actin ([l9l ). 



When discussing the effect of force on a single actin or microtubule fil- 
ament, one important issue is the buckling of the filament. Since actin 
filaments have much smaller persistence length Ip than microtubules, actin 
filaments buckle easily under external force. Our approach is appropriate to 
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describe experiments like that of ref (j34l ). where very short actin filaments 
are used. The length of the filaments must be smaller than the critical length 
for buckling under a force /. This length can be estimated as = -i\\J n/ f 
with a hinged boundary condition, where k = IpksT. Taking Ip ~ 17/im, 
we estimate 1}^ ^ 828nm at / = IpN. Our discussion of the force will be 
applicable only for filaments shorter than 



Collapse time 

In this new section, we shall study experimentally relevant questions such 
as the mean time required for the ATP cap to disappear or the mean time 
required for the whole filament (ATP cap and ADP subunits) to collapse to 
zero length. We are interested in the conditions for which these times are 
finite. Below we address these questions. 



Cap collapse in phases I or II 

The dynamics of the cap corresponds to that of a ID biased random walker 
with a growth rate U and a decay rate Wt + R- Here, we calculate the 
mean time required for a cap of initial length kd to reach zero length for 
the first time. We assume that there is a bias towards the origin so that 
Wt + R > U . This time is nothing but the mean first passage time for 
the biased random walker to reach A; = 0, starting from an arbitrary site 
k in phases I or II. According to the literature on first passage times, the 
equation for is (0, 0, Eil) : 



Un+i + {Wt + R)Tk-i -{U + Wt + R)Tk + 1 = 0. (39) 

When Wt + R > U, this recursion relation can be solved (see Appendix B) 
with the condition that Tq = 0, and we obtain 

^ (40) 



Wt + R-U 



This corresponds to the time the random walker takes to travel a distance k 
at a constant velocity — J. Note that the mean first passage time becomes 
infinite in the unbiased case when J = or if the bias is not towards the 
origin i.e. when Wt+R < U (which would correspond to an initial condition 
in phase III) 

One can also define an average of the mean first passage time with respect 
to the initial conditions. Averaging over k and using equation (jlSh . one 
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obtains 



{Tk) 



U 



(41) 



{Wt + R-UY' 



The same time can be recovered by considering the average time associated 
with the fluctuation of the cap: 



This time may be related to the catastrophe rate in the foUowing way. In 
ref. 

(0) 

, the catastrophe rate is defined as the total number of catastrophes 
observed in an experiment divided by the total time spent in the growing 
phase. Since the growing phase ends when the cap disappears for the first 
time, we interpret similarly l/{Tk), as an average collapse frequency of the 
cap. 

Filament collapse in phase I 

Now we consider the dynamics of the filament length which is described 
similarly by a 2D biased random walk converging towards the origin. Here 
we investigate the mean time T„_fc required for a filament with an initial 
state of n ADP subunits and k ATP subunits to reach zero length for the 
first time with an initial condition inside phase I. Again, this is the mean 
first passage time now in a 2D domain (in the n — k plane, as shown in 
Fig. [1]) to reach the origin (n = 0,k = 0) starting from an arbitrary n and k. 
This mean first passage time r„ ^ obeys the following set of equations (jsoj). 
When k > 0, for all n, the equation is 



UTn,k+i + WrT^^k-i + RTn+i^k^i -{U + Wt + R)T^,k + 1 = 0. (43) 



and we also have the condition Too = 0. 

The simplest way to solve these equations is to guess by analogy with 
the ID case that the solution must be a linear function of n and k. This 
leads to a simple ansatz of the form Tn^k = An + Bk, which in fact gives 
the exact result as can be shown rigorously. Substituting this in Eq. |33] and 
Eq. H¥l we can solve for unknowns A and B. This leads to 



{Tk) 



(42) 



2D^ 



For k = Q and n > we have a special equation 



UTn^l + WoTn-lfi -{U + WD)Tnfi + 1 = 0, 



(44) 




(45) 
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Figure 7: Mean time taken by a filament of initial length {n + k)d to collapse 
to zero length. Curves are given by eq. US] and points are obtained from a 
Monte Carlo simulation for different values of n and k. From bottom to top 
(n, k) = (990, 10), (750, 250), (500, 500) and (200, 800). 

where vu is the velocity of the intermediate phase given by Eq. [25j Note 
that vn < here since the initial condition is within phase I. 

We first examine some simple particular cases of Eq. [JSj As we approach 
the intermediate phase boundary vn 0, Tn^k ^ oo as expected. When 
— oo, Tji k — k/i^Vx -\- R — — which is the cap collapse time 
calculated in the ID case. When Wd = oo, the whole filament collapses 
immediately after the cap has disappeared for the first time i.e. after a time 
Tfc. When R ^ oo, ATP subunits instantaneously become ADP subunits 
and we obtain another simple result Tn^k = {n + k)/{WD — U). We have 
also compared the prediction of Eq. US] with Monte Carlo simulations in 
Fig. [7] and we have found an excellent agreement. 

We can also define an average of the above mean first passage time where 
the average is performed over initial lengths of cap and un-hydrolyzed region. 
Averaging over k and n in Eq. [35] we obtain 



The inverse, l/(r„ fc), can be called the collapse frequency of the filament. 
The filament collapse frequency and the cap collapse frequency are shown in 
Fig. [8] as a function of U and C for the cases of actin and microtubule using 



{Tn,k) 



U{R^ + WdR + WtR + Wl) 



(46) 



{UR - WdWt - WdR + UWd) 
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C (|xM) C (nM) 

0.095 0.112 0.129 8.125 8.313 8.5 




Figure 8: Collapse frequencies for actin and microtubule(MT) as function of 
the rate of addition of monomers U (lower x-axis) and of the concentration C 
(upper X-axis): Solid line represents the collapse frequency of the filament 
given by, l/l^n,k)i the inverse of the time given in Eq. H6l Dashed line 
represents the collapse frequency of the cap given by, 1/ (T^), the inverse of 
the time given in Eq. W\\ 

parameters of Table [H Both frequencies are close to each other because 
the rate Wd is large compared to other rates (see Table [I]) . This figure also 
shows that as the frequency of collapse is increased, the rate U decreases and 
so the filament length is decreasing, which is expected 

0). 

The behavior 

of the collapse frequency as function of the growing velocity in the absence 
of force agrees with (26. ). The decrease of the rate of monomer addition 
is in practice caused by either the application of a force or a lowering of 
the concentration. Thus the application of force may be seen as a general 
mechanism to regulate the dynamic instability. 

In the opposite limit, when R/Wd — > one can also understand the 
result physically from the following argument. When the filament collapses, 
the first event is the disappearance of the cap and therefore the first contribu- 
tion to the collapse time is the mean time required for the cap to disappear, 
Tfc, as obtained from Eq. (|40p . Once the cap has disappeared, assuming that 
Wd is very large, ADP subunits start depolymerizing until the next ATP 
subunit addition takes place. The mean time needed for an ATP subunit 
addition to take place is 1/U . Once an ATP subunit is added, one has to 
wait an average time of Ti for the cap to disappear again. This cycle of 
ATP subunit addition and depolymerization repeats many times. The num- 
ber of times this cycle occurs, starting with a filament of n ADP subunits, 
is roughly n/{^jf-). But one also has to take into account the increase in 
ADP subunits as a result of ATP hydrolysis, which is done by subtracting 
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RTi from n/(^^). This leads to the fohowing approximate expression for 



' n.k 



nd 



n 



( Wp 
\ u 



RTi) 



-VII 



+ 



kd_ 



(47) 
(48) 



where J = {U — Wt + R)d is the cap velocity in the rapidly growing phase. 
This solves Eq. |43]and Eq. [44]in the limit R/Wd — > and agrees reasonably 
well with the Monte Carlo simulations. 



Conclusions 

In this paper, we have studied a model for the dynamics of growth and 
shrinkage of single actin/microtubule filaments, taking into account the 
ATP/GTP hydrolysis which occurs in the polymerized filament. We find 
three dynamical phases with different properties of the ATP/GTP cap and 
the filament: a bound phase, an intermediate phase and a rapidly growing 
phase. For each phase, we have calculated the steady-state properties of the 
non-hydrolyzed cap and of the filament and we have investigated the role of 
an external force (f) applied on the filament during polymerization and of 
the monomer concentration (C), leading to a f-C phase diagram. We have 
also calculated the collapse time which is the time needed for the cap or the 
filament to completely depolymerize. 

In the bound phase, both the size of the filament and the size of the 
cap are finite. It must be noted that in batch experiments, where the total 
(free+bound) concentration of monomers is a constant, the steady state 
is always the bound phase. If the total monomer concentration is large 
enough, the filament length is very large and the concentration in the bound 
phase will be just below the critical concentration. The filaments exhibit 
large length fluctuations in the bound phase corresponding to the dynamic 
instability observed with microtubules. An important conclusion is the role 
of ATP/GTP hydrolysis in determining the stall force of the fllament: the 
effect of hydrolysis always reduces the value of the stall force. Dynamic 
instabilities are in general not observed with actin filaments when no force 



is applied to the filament. However, in recent experiments (|2ll . l22l ). it 
has been found that the presence of ADF/Cofilin leads to a bound phase 
even at zero force where actin filaments exhibit large length fluctuations. 
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In the absence of binding proteins, a natural way to regulate the dynamic 
instability and the length of the filament is through the application of force. 

The intermediate phase is a phase where the filaments grow at a constant 
velocity with a finite ATP/GTP cap. This is the phase which is in general 
observed in a cell. There are large length fluctuations of the cap in this phase 
but the length fluctuations are not as large as the average length. Thus, there 
is no true dynamic instability in this phase but we give predictions for the 
typical time of the collapse of the cap. Finally, the growth phase corresponds 
to the case where both the filament and the cap grow at a constant velocity. 

One of the limitations of our work is that we considered a single protofil- 
ament. This does not seem to be an important issue for actin filaments 
where the length difference between the two protofilaments is always small 
and of the order of an actin monomer. For microtubules, the detailed poly- 
merization mechanism seems to be very complex and this certainly plays an 
important role on the way the force is distributed between protofilaments 
(j42l ). Recent experiments have also considered the maximum force that can 
be generated by a bundle of parallel actin filaments (j34l ) . Our results raise 
interesting questions about the role of ATP hydrolysis in this case. 

Our work could be extended in several directions. In the biological con- 
text, actin or microtubule polymerization is regulated by capping proteins 
and it would be important to understand quantitatively the regulation mech- 
anism and to incorporate them in our model. This is for example the case 
for motors of the Kin-13 family that have been found to interact with micro- 
tubules and induce filament depolymerization 043). So far, we have mainly 
considered the polymerization kinetics, but in general, there is a complex 
interplay between the mechanical properties of the filaments and the poly- 
merization kinetics, which we plan to explore in future work. 
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Appendix 

A. Calculations using the generating function approach 

Summing over n and k in Eq. [1] and using equations [21 [3] and [71 one obtains 



dG{x,y,t) 
dt 



U{x-1) + Wt 



Wt 



e-0 



+R{--1 

y 



G{x,y,t) 



1] +R[- -1] +Wd{1-- 



+ WD{i-^]p{o,o,t). 



Foix,t) 

(49) 



Similarly one can also write down equations for and Hn- This equation 
contains Fq, which is coupled to all the For A: > 0, 



dFk{x,t) 



dt 

and for A; = 0, 

dFo{x,t) 
dt 



UFk^iix, t) + {Wt + Rx) Fk+iix, t)-{U + WT + R) Fk{x, t), 

(50) 



Wd{^-1] - U] Fo{x,t) + {WT + Rx)Fi{x,t). (51) 



Solving this set of equations we shall derive a formula for G{x,y,t). From 
G{x,y,t), we calculate the following quantities: 
The average length 



(/) = [(n) + (k)] d = d 



the velocity of the filament 



dG{x,l,t) 
dx 



d 



x=l 



dG{l,y,t) 
dy 



y=l 



d{l) d f dG{x,x,t) 
V = lim — — = d lim — ; 

t^oo dt t^oo ox \ dt 

and the diffusion coefficient of the filament length 



(52) 



(53) 



x=l 



d^ lim 



ld^fdG{x,x,t)\ , 1 d fdG{x,x,t) 
2d^ 



+ 



dt J ' 2dx 
dG{x,x,t) \ d_ f dG{x,x,t) 
dx I dx \ dt 



dt 



(54) 



;r=l 
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The average velocity of the cap is 



J = ,n.njf> = lio.|,(-'g(',f'>) , (55) 



t^oo dt t^oo dy \ dt y 
and the diffusion coefficient of the cap is 

t^oo2dt ' ^ ' ' 



cP hm 

t— ►oo 



1 92 fdG{l,y,t)\ 1 d (dG{l,y,t 



2 5yH J 2 9y V dt 



dG{l,y,t)\ d /dG{l,y,t 



dy J dy \ dt 



(56) 



y=l 



Calculation of F{x = l,t ^ cxd) in phases I and II 



In the steady-state (t 00), the cap distribution in phases I and II becomes 
time-independent and hence [dFk/ dt)x=i = 0. In this case Eqs. [50l and [5T] 
can be written, for A; > 0, as 

= UFk-i + {Wt + R) Fk+i -{U + WT + R)Fk, (57) 
and for A; = 0, 

= {Wt + R)Fi- UFq. (58) 

where we denote for short, F^ = Fk{x = l,i — > 00). The solution of Eq. [ST] 
is of the form F^ = q^Fo. If we substitute this back into Eq. [571 we get a 
quadratic equation in q 

{Wt + R)q^ -{U + WT + R)q + U = 0. (59) 

The two solutions are q = U/{Wt + R) and q = 1, but we can rule out 
q = 1 using the normalization condition YlT=o^><^ ~ ^- phases I and II, 
Wt + R > U and therefore q < 1. Using the normalization condition, we 
obtain 

Fk = {l- q)q\ (60) 

which is Eq. [T3l 

Calculation of G(x, y) in the bound phase (phase I) . 

We now explain how to calculate G(x, y) in the bound phase, using a tech- 



nique of canceling apparent poles (j4Cl . l4ll ) . Since we are interested in the 
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steady-state properties of the bound phase, the time derivative of G on the 
left hand side of Eq. |39] is zero, which leads to 

_ Fo(x) [Rxjy -x) + WTxjy - 1) - Wpyjx - 1)] - WpPjO, 0)y{l - 
x[-Uy^ + iU + WT + R)y-Rx-WT] 

(61) 

where Fk{x) and P(0, 0) are unknowns. By definition 
G(x,y) = ^ J]P(n,A;)xV, 

n>0 fc>0 

since the P(n, k) are bounded numbers, G{x, y) is an analytic function for 
< |x| < 1 and < |y| < 1. To guarantee the analyticity of the function 
G{x, y), the zero of the denominator of Eq. [HI 

y = y- = :^(u + WT + R- ViU + Wt + Rf - W{Wt + Rx)) , 



2U 

must also be a zero of the numerator. This implies that 

p . X WDP{0,0)y.{l-x) 

' Rx{y^ -x) + WTx{y. - 1) - Woy-i^ - 1) ' 

The normalization condition, namely G(x=l,y=l)=l, then fixes the value 
of P(0, 0) as 

^ ' ' Wd \Wt + rJ 

After substituting Eqs. [62] and [63] into Eq. [611 we obtain the expression of 
G{x, y) given in Eq. [TO] 

Velocity and diffusion coefficient in the Intermediate phase (phase 
11) 

We recall the definition of Ffc(x, t) given in Eq. [7] 

Ffc(x,t) = J^P(n,A:,f)x", (64) 

n>0 

and we recall that -F^ with no argument is a short notation for Fk{x = l,t 
oo). From this, we introduce 



ak{t) 
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so that < n{t) >= Efc>o -Pl^' = Efc>o"fcW- 

By taking a derivative with respect to x in Eqs. [50land[5H one obtains 
the equations of evolution of ak{t): for A; > 

= Uak_iit) + {WT + R)ak+iit)-iU + WT + R)ak{t) + RFk+,{x = l,t) 

(66) 

and for A; = 

= -WdFo{x = l,t)-Uao{t) + {WT + R)ai{t) + RFi{x = l,t). (67) 

at 

As shown in (jj^ ). there is a solution of these recursion relations in the long 
time limit in the form of ak{t) = M^t + B^, where Mk and are time 
independent coefficients. After substituting this equation into Eqs. [6611671 
and separating terms which are time-dependent from terms which are not 
time-dependent, one obtains separate recursion relations for and Bi^.. 
The recursion relation of M^, is identical to that of obtained in Eq. [57H58[ 
Using Eq.[60l the solution can be written as Mj. = vii{\ — q)q^ / d = vjjFk/d. 
The recursion relation of Bk is for A; > 0, 

Mk = UBk-i + {Wt + R)Bk+i -{U + Wt + R)Bk + RFk+i, (68) 
and for A; = 0, 

Mo = -WdFo-UBq + {Wt + R)Bi + RFi. (69) 
These recursion relations can also be solved with the result 



Bk = B^q^ + 



Wd{1 



Wt + R 



kq''. (70) 



To characterize the intermediate phase, it is convenient to rewrite the 
evolution equation for the generating function G{x, y, t) of Eq. Housing the 
fact that i-*(0, 0, t — > oo) = in this phase, in the form of an evolution 
equation for G(x, t) = G{x, x, t) as 

^^^ = a(x)G(x,t)+/3(x,t), (71) 

where a{x) = U{x - 1) + Wt{1/x - 1) and I3{x,t) = (1 - 1/x){Wt - 
WD)FQ{x^t). With this notation, the velocity defined in Eq. [53l is 
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where the prime denotes derivatives with respect to x. Substituting the 
expressions of a{x) and f3{x^ t) into this equation, it is straightforward to 
obtain the velocity v = vu characteristic of the intermediate phase which is 
Eq. [25l Similarly, using Eq. [5l] and Eq. [Til the diffusion coefhcient can be 
written as 



D = 
where 



— lim 

2 t^oo 



a"(l) +/?"(!, t)+a'(l) + /?'(l,t)-2G"(l,t)/?'(l,t) , (73) 



G'iht) 



dGix,t) 



dx 



vt + 



1 



x=l 



--{n{t)) + {k), 
U + Wpq 

Wt + R _ 



(74) 



After substituting Eq. [74] into Eq. [73] and simplifying, the terms linear in 
time and the term containing the unknown parameter Bq cancel out in the 
expression of the diffusion coefficient, and we finally obtain D = Djj which 
is given in Eq. | 



B. Cap collapse time 

In order to solve the recursion relation for of Eq. 
transformation defined by 



T(z) 



k>0 



After using the initial condition Tq = 0, we obtain 

[l + UTi{l- z))z 



T{z) 



U{l-zY{q-^ -z) 



we perform a z 
(75) 



(76) 



where = {Wt + R)/U . By definition, T{z) is analytic for all values of 
\z\ > 1. Since we are interested in the case Wt + R > U, the numerator 
in Eq. [76] must vanish to ensure that T{z) is analytic at z = q~^. This 
condition determines the unknown Ti = 1/{Wt + R — U). Now can be 
obtained by an inverse Z transform as 



1 
27d 



T{z)z^-^dz 



1 
27ri 



z dz 



k 



C/(l-z)2(g-i -1) Wt + R-U' 



(77) 
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